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ABSTRACT 

We consider the application of logarithmic conformal field theory in finding solutions to 
the turbulent phases of 2-dimensional models of magnetohydrodynamics. These arise 
upon dimensional reduction of standard (infinite conductivity) 3-dimensional magne- 
tohydrodynamics, after taking various simplifying limits. We show that solutions of 
the corresponding Hopf equations and higher order integrals of motion can be found 
within the solutions of ordinary turbulence proposed by Flohr, based on the tensor 
product of the logarithmic extension cq^ of the non- unitary minimal model c^x. This 
possibility arises because of the existence of a continuous hidden symmetry present in 
the latter models, and the fact that there appear several distinct dimension -1 and -2 
primary fields. 
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1 Introduction to Polyakov turbulence and Flohr 
solution 



The word turbulence is used to describe the irregular and fluctuating motion that is 
seen in many cases of gas and liquid flow. We can have a qualitative description of the 
flow by a dimensionless number, the Reynolds number Re, which measures the relative 
strength of inertial and viscous forces acting on the fluid. Turbulence is associated 
with flows having very high Reynolds number. Formally we can take Re — > oo and this 
is the case of the fully developed turbulent flow. The Navier-Stokes equation describes 
these flows in 2 and 3 dimensions. One of the main differences between inviscid flows in 
each of these dimensions, is the existence of an infinite number of conserved quantities 
in 2-dimensions, and it is these and their corresponding conserved fluxes that make 
two dimensional turbulence particularly interesting. Indeed there has been a growing 
body of research on the subject since the classic work of Kraichnan [1]. 

A fresh approach to the problem of understanding turbulent flows in d — 2 was 
presented by Polyakov [2,3], who, motivated by the success of Conformal Field Theory 
(CFT) in describing 2D critical phenomena, suggested that the methods of CFT can 
help in this task . The central point in the analysis was the Hopf equation, written 
here in terms of the vorticity u = e a pd a vp, a, (3 = 1, 2 



p=i 

where the time derivative of uj appearing in eq(|I|) is given in terms of the Navier -Stokes 
equations 



where, assuming incompressibility, the velocity is expressed in terms of the stream 
function ip as v a = e a pdpip 

The Hopf equations express n-point function in terms of the n + 1 point functions. 
They can't be solved directly and the only way so far was to take ad-hoc approximations 
and try to find solutions. In contrast, with the Conformal approach we try to solve 
the Hopf equation exactly by the ansatz that the flow has conformal symmetry in the 
situation where it is fully developed. 

The assumption is that there exists a range of scales such that the affect of the 
viscocity term is very small - the so called inertial range. So we insert u = —d 2 ip in 



n 




(1) 



d t uj + e a pd a ipdpuj = vd a d a uj 



(2) 
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the Hopf equation, postulating that ip is a primary operator of the CFT we are looking 
for. Because of this, one is naturally lead to introduce point splitting regularization 
in order to evaluate products of fields defined at the same point. Using the operator 
product expansion (OPE) for ip with itself 

M x M = \<t>] + •••• (3) 

where is the conformal field with the minimal dimension appearing in the OPE, and 
+.... indicates next to leading order terms, the Navier-Stokes equations imply that 

Co = -|a| 2 ( A ^ 2A ^£0, C = [L_ 2 L\ - L_ 2 L\] (4) 

appears inside the correlators in the Hopf equation, where a is the u.v. cutoff associated 
with point splitting. In principle, a solution to the Hopf equations can be found if the 
operator on the right hand side in eq(f|) is vanishing. However, this requirement by 
itself is not restrictive enough to allow us to find explicit solutions. To place further 
restrictions, one has to take into account the spectrum associated with turbulence. 
One of the predictions of the conformal approach is that the fc-space kinetic energy 
density E(k) (defined by / dkE(k) = (v a (k)v a (—k) ) ) implies a power law spectrum 
with E ~ /c 4A * +1 . In the classical turbulence studies, the value of the exponent of k in 
the above relation was one of the first things to be studied. There are different opinions 
for its value but the work of Kraichnan [1] indicated that in the case of the inverse 
cascade scenario, (similar to the Kolmogorov one in d — 3), E ~ k~ 3 which is based on 
the constancy of the enstropy flux. So it is natural to postulate this cascade scenario 
in order to increase our restrictions on the possible conformal field theory solutions. In 
this case we have 

(u(z,z)u(0) ) ~ const A + A^, + 3 = O (5) 

Trying to find solutions of these constraints, Polyakov considered the possibility 
2A^, < A^, which 'trivially' satisfies the Hopf equations in |a| — > limit, and found 
solutions within the non-unitary minimal models. In fact hundreds of solutions can 
be found in this way [4-9], which is perhaps an indication that this approach is not 
encapsulating all the physical aspects of turbulence we would like. Worse still, Falkovich 
and Hanany [10], showed that a solution that has the Kraichnan spectrum can not be 
found within the minimal models. Because of the constancy of all the higher order 



2 



fluxes, they found that 2A^ = A^ = —2 and this case is not found in the minimal 
models. Also in the minimal model case the actual prediction is that all fluxes vanish 
so the situation is made more difficult. Such difficulties indicate perhaps that the 'usual' 
Conformal approach might describe certain kinds of flows, but not the turbulent ones. 

In an attempt to directly deal with these issues, Flohr [11] has suggested a solu- 
tion that does not suffer from these problems and is based on Logarithmic Conformal 
Theories (LCFT) [12 -16]. This ability to avoid the aforementioned difficulties and yet 
still remain within the conformal approach is clearly due to the novel characteristics 
of LCFT's. 

One such important characteristic, is the existence of pairs of primary fields C(x), D{x) 
in the spectrum having the same dimension and with respect to which the Virasoro 
operator L is non-diagonal i.e. L \C) = Ac\C), L \D) = Ajj\D) + |C). 

The two point functions of the logarithmic pairs C(x), D(x) have the following form 
(see e.g. [13] ) 



(C(x)D(y)) 



d 



(x-y) 



2A r 



(D(x)D(y)) = d -—(-2e\n{x-y)+d) 

[x-y) 

(D(x)D(y)) = (6) 

where d,e are constants. As pointed out in [13], the first correlation in eq(f|) together 
with the invariance of the correlator under exchange of C(x) and D(x) implies a strong 
constraint on the allowed values of the dimensions Ac, namely that it be integer.^] 

This implies that the logarithmic operators C (x) (or rather their chiral components) 
are associated with a hidden continuous symmetry, although unlike the more usual 
current algebras, the Schwinger term is vanishing. 

This hidden continuous symmetry associated with LCFT's is an important ingre- 
dient in Flohr's solution of the Hopf equations ([1]). 

In searching for possible non-unitary LCFT solutions to turbulence, (the non- 
unitarity is a consequence of the flux constraint (H)), the logarithmic extensions of 
the minimal model CFT's is perhaps the first place to start since a certain amount 



3 See also [26] for additional clarification of this point 
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of knowledge already exists about their properties [16]. In fact it is the c Pj i models 
that have been studied most. As we have seen above, the Kraichnan spectrum implies 
that there should be A = — 1 and —2 states in the spectrum of our LCFT. This is not 
possible in a single copy of c p> \ but as shown by Flohr [11], since c 6j i does contain a 
A = — 1 field, the tensor product c 6j i®c 6) i can contain both. The additional advantage 
of the tensor product model is that it naturally contains Polyakov's Bose condensate 
[3] within the solution. 

C6,i has 17 primary fields, 5 of them are "logarithmic". The dimensions of the fields 
are A 1 >r = y-((p — r) 2 — (p — l) 2 ) with p = 6. 

Because of the tensor product there is a degeneracy in the way we can define the 
states with dimension —1,-2 . In fact this will prove to be a nice feature of the 
solution, which we shall exploit later on when considering the MHD case. 

We define the A = — 1 and A = —2 fields, ipi,<pj' as follows (with $(aia') = 
(<3>a <8> $A')(z), and <&a(#) a primary field in c 6 ,i ) 

V>1 = ^(^(-l|-0) + ^(-0|-l))> ^2 = ^($(_i|_o) + $(_0|-l)) 

0i = <t>2 = -(%ai-i) + %i h i)), 03 = $(_i|-i) (7) 

It is easy to see that x ipj = J2k' C^^fc' + ••• >*> 3 — 1)2 and fc' = 1, 2, 3, where the 
OPE coefficients C^i can be obtained from the fusion rules of cq^ LCFT given in [16]. 
We expect the physical fields to be given by linear combinations of the fields ipi and 

4>f, i.e. ijjp h = Ei ctiifri and <j) ph = Ej' Pj'<f>j> with 

4> P h x Vpft = P h + (8) 

In order that the physical fields ip p h satisfy eq(|8|), the coefficients ai and /3j> must 
satisfy ^(y-jCijk 1 = ftw which essentially determines the coefficients f3k> once the ai 
are known. 

From the general form of ipph one can determine that its fc-space two-point function 
has the structure 

(MQM-k)) = + ^s(k) + .... (9) 

with the coefficients ao ~ (I) and ai ~ (4>ph)- A notable feature as pointed out in [11] , 
is the presence of the so called Polyakov condensate term (proportional to S(k)) which 
is naturally contained within the solution. 
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Flohr considered a more general expansion of the physical fields i/j p h and (j) p h involv- 
ing non-diagonal combinations of left and right-moving states. This leads to non-zero 
spin contributions to ip p h and 4> p h that may naturally encode the effect of stirring forces 
in the fluid. The general expansions are of the form 

4> P h(z,z) = ^ph(z) x ip ph (z) + Y U(z)xV(z) 

Au+A v =-2, A v +A v eZ 

ph (z,z) = <j> ph {z) x (t) ph {z) + Y U'{z)xV'{z) (10) 

A c/ /+A v ;=-4, A u ,+A v ,eZ 

where the integer dimension fields U(z), V{z), U'(z), V'(z) are those occuring in C6,i <8> 
c 6) i. This leads to the following expression for ip ph 

ip ph (z,z) = ai$ ( _i| )(z)®$(_i|o)(^) + a 2 $(o|-i)(^)<2)^(o|-i)(^) 

+ a3$(_i|6)(*) ® $(-i[o)(*) + a 4$ ( o|-i)(^) ® $(6|-i)(*) 

+ a 5 $(_i|_i) (-z) (8) $(o[o)(«) + O6^(o[o)(-z) ® $(-i|-i)(^) 

+ a7$(o|o)(*)®$(6|o)(s) + '" ( n ) 

There is a similar expansion of (f> p h{z, z). The non-diagonal terms in eqflTTD induce 
further terms in the two-point function eq([|) that are proportional to the expectation 
values of chiral fields, and involve the derivative of fc-space delata functions [11]. A 
crucial aspect of the model is the existance of a hidden continuous symmetry [13] 
associated with the integer dimension field \l/ Ai 5 • I n fact this guarantees that it is a 
solution to the Hopf equations. The operator Q(z) = L X [L_ 2 -^ 2 -i — L _ 2 L 2 _i\(i> p h(z) 
generates a continuous symmetry of the LCFT and satisfies the condition [] 

(fi(z)$ 1 ( 2l )$ 2 (^)....$ n (z„)) = (12) 

In ending this section, it is perhaps worth mentioning the reason one finds non- 
unitary CFT in connection with turbulence, may in some sense be expected, since 
the turbulent state is not an equilibrium state. Even if the conformal solutions don't 

4 As argued in [11], The presence of L L in f2(z) is needed to ensure that only the ^>a 1[S part of 
4> p h enters the Hopf equations. This follows from certain properties of correlation functions of LCFT 
discussed in [15] 
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describe the turbulent state exactly, we might at least expect them to describe flows 
out of equilibrium. 

Motivated by the new possibilities that LCFT has brought to conformal turbulence, 
in the next section we want to extend Flohr's construction and look for solutions of the 
generalized turbulent MHD systems in d = 2 that arise naturally when dimensionally 
reducing standard ideal MHD in d = 3. 



2 Logarithmic solutions to MHD turbulence 

The long wavelength limit of 3-dimensional plasmas are well approximated by the 
equations of Magnetohydrodynamics (MHD), which in the limit of low frequencies, 
low temperature and incompressibility take the simplified form [17] 

V-v = 
V-B = 

9 f v + vVv = — -VP+— JxB (13) 

Pm CPm 

V x (v x B) = d t B 

V x B = —J 

c 

Here v, B, J are the velocity field, magnetic field and electric current respectively P 
and p m are the pressure and mass density. After simple manipulations, we can recast 
these equations into a more symmetric form 

V-v = 
V-B = 

V x (d t v + v • Vv) = — V x (B ■ VB) (14) 

A-np m 

(B • V)v — (v ■ V)B = d t B 

In order to obtain an effectively two-dimensional theory from eq(|I~4"D, one simple 
proposal is to dimensionally reduce the model by requiring that all fields be independent 
of the z-coordinate, i.e. d$v = and c^B = 0. It follows that the first two components 
of v and B can be written as 

B a = e a /3dpA, v a = e a pdpip, a, (3 = 1,2. (15) 
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One can define, in analogy with vorticity u, the "magnetic vorticity" 



ft = e aP d a Bp = -d a d a A, (16) 
and the two-dimensional operator 

A = e aP d p Ad a (17) 

after which the 2-D MHD equations take the form [18]: 

Co + £ a 0df3tpd a u) = — ^ — Afl 
A + tapdpipdaA = 

V + e aP d^d a V = — ^ — AB (18) 

4vrp m 

B + e a pd^d a B = AV 

Solutions to eqs(^) have been studied in the context of minimal model CFT, in 
two simplifying cases, namely parallel flow conditions where B = V = [18], and 
perpendicular flow for which A = V = [19]. In both cases two of the four equations 
in (p~8| ) are automatically satisfied and the remaining equations to be studied are: 

^ + tattdalpdaU = Afl 

4?rp m 

A+EapdfrfdaA = (19) 

for parallel flow and 



Co + eapdpipdaUJ = 

B + e aP dpipd a B = (20) 

for perpendicular flow. 

We want to discuss the behaviour of various fluxes when one considers dissapative 
forces. In the case of ordinary fluids this means reintroducing the viscosity v. In the 
MHD case there is an additional dissapation parameter namely the conductivity o. 
Keeping v = for the present but a finite, we find for the equations of motion of 
perpendicular flow [19] 
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uj + € a pdpil)d a uj = -— (-^—d a d a B - -[e aP d^d a B\) (21) 

p m Ana c 



c 2 



B + e af3 df3^d a B = —d a d a B 

Finite viscosity v would add the standard dissapative term to the rhs of the ui 
equation of motion but leave that of B unaffected. Hence viscosity and conductivity 
are responsible for the decay of higher point enstrophy and magnetic fluxes, H n (k,t) 
and B n (k,t) where 

J dkH n {k,t) = J d 2 x(u n {x)), J dkB n {k,t) = J d 2 x(B n {x)) (22) 

For infinite conductivity but non-vanishing u, B n is time independent (so 5-fluxes 
are conserved) but H n satisfies 

ra- 1 

H n ~v k 2 h n (k,t)dk (23) 
Jl-1 

where we have introduced the I.R. cutoff L. Dimensional analysis implies that the 
viscosity v and ultraviolet cutoff a are related via v ~ a 2A *~ 2A '*. A further analysis on 
the viscosity dependence of H n {k) implies [5] 

Hn ~ Z/ [(n_1)(A ^ +1)+A '* +2l/(A ''>~ A ^ ) (24) 

Demanding that the higher order fluxes H n have the same scaling properties in- 
dependent of n determines the scaling dimensions of both and ip to be = 
— 1, = —2 and immediately gives rise to the Kraichnan form for the energy spectrum 
E(k) ~ k~ 3 [5]. 

In the present case of MHD we can also imagine considering the situation where we 
have inviscid flow but with a finite. In this case we can perform a similar dimensional 
analysis as before. The relation between a and a can be determined from the B 
equation of motion (20) to be a^ 1 ~ a (2A B -2A x ,-2)^ w h ere ^' j s defined to be the 
minimal dimension field in the OPE B x \&. The fluxes H n and B n then satisfy 

H n ~ ^a- 1 [d 2 x(uj n - 1 d 2 B} + ... 
B ~ a' 1 J d 2 x(B n ' 1 d 2 B) (25) 



8 



It follows that the operators u and B have the same dimension and hence A x > = 
A^ = —2 for consistency. The a dependence of B n can be obtained as 

B n ~ (r Kn-2)A B +A x ,+2]/(A B -A x ,) ^ 

Hence it follows that the scaling behaviour of B n is independent of n if A# = 
and A x i = —2, and we see the condition on A x > agrees with the one found above by 
consistency. 

In studying the possible solutions of the inviscid ideal MHD equations, in the per- 
pendicular flow regime eq(^), we have to identify the dimension zero field B with 
certain primary operators. We want to look for solutions of the resulting Hopf equa- 
tions similar to the logarithmic conformal field theory solutions considered by Flohr in 
the case of normal fluids. An obvious motivation for this is that the LCFT cq^i ® 5q ; i 
considered in [11] solves the Hopf equation involving u, which is one of the equations 
in (^). The hope is (and we shall see that this is indeed possible) that one can identify 
the B field with primary fields within this particular LCFT in such a way as the Hopf 
equations associated with the second equation in ( |2Tip can be satisfied. Even though 
this is in the context of perpendicular flow, one should expect that similar ideas would 
apply to the parallel case as well. 

Indeed, considering the equation of parallel flow, eqflHf), one obtains after perform- 
ing the regularization and OPE's, 

<j = \a\ 2 ^- 2 ^&p - |a| 2 ^- 2A ^£0 = 

A = |a| 2{A *- A ^ A ^£x = (27) 

where in eq(|27[), C = [L^L 2 ^ — L^L 2 ^} and A x A = [p\ H ,ip xip = [(f)] H and 

A x if) — \x] + ■ • ■■ Here (p, <f> and \ are minimal dimension fields. 
Perpendicular flow corresponds to the equations 

u) = -H 2 ^- 2A ^£0 = O 

B = -\a\ 2{A x'~ A ^ AB) C X ' = (28) 

where \' is the minimal dimension field in the OPE B x if) = [x'] + • • •. 

Now we have seen previously that flux conditions dictate that A# = 0. Also 
one finds that on dimensional grounds A^ = 0. Thus both these fields have to be 
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constructed from A = fields in the c 6 ,i x c 6j i LCFT. Just as in the case of the 
A = — 1, —2 states, there are several possible primary fields of dimension zero, which 
we call by, i' = 1, 2, 3, where 



h = $ ( 6|6), h = -($(i|_i) + $(-in)) 

h = (29) 

(There are other possibilities involving tensor products containing the identity op- 
erator / which we do not include as they are in some sense trivial). 

As in the case of ordinary turbulence, we expect that the physical components of the 
magnetic field B and A will be given by some linear combination of the bi fields. Under 
these circumstances, in verifying either of the two sets of Hopf equations following from 
eqs(^j], ^S]) we need to calculate the OPE's by x by and by x ipj. In order to obtain these 
expansions we need some further fusion rules of C6,i beyond those considered in [11]. 
Specifically we need those involving the primary field [1] with the set [0], [—1], [—1], [1], 
which can be obtained from the general rules presented in reference [15], 

[1] x [0] = [-|] + [-1] + [-1] + [-|] + [6] + [1] + I 7 -] + [4] 

[l]x[-l] = [-|] + [0] + [1] + [Z] + [4] 

[l]x[-l] = [-1] + [-| + [0] + [1] + [I] + [4] 

[l]x[l] = [0] + [-|] + [-1] + [-|] + [1] + [I] + [4] (30) 

where the right hand side of the OPE's in eq(|30D is to be understood modulo multiplic- 
ities and subject to the linear relation amongst certain conformal families as explained 
in detail in [16]. Using these relations and those defined earlier in section 1, one can 
obtain 



bx x *0i = 402 

bx x tp 2 = 20i 

b 2 X Ipx = -01 + 02 + -03 

b 2 X 1p 2 = 02 + 01 
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1 
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-d 
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)i -1 — a 
1 






03 x W2 — 
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(p2 


1 , 
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2 
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01 










1 




1 




bi x b 2 = 
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>1 + 02 
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03 


61 X 63 = 


01 


+ 02 
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1 






b 2 xb 3 = 


-c 




>2 




4 
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1 






b 3 xb 3 = 


-4 




>3 + 


^02 


2 









An important feature of the OPE in eq(|3TD is that the leading terms on the right 
hand side (i.e. those with smallest dimensions) are those of [—1] and [—1], except in 
the OPE [1] x [—1] where the lowest dimension is —2/3. Despite this exception one 
finds that the leading order representations on the rhs of eq(|30|) are just given by linear 
combinations of the 0j fields defined earlier. 

In the case of parallel flow the solution includes a field with dimension —1 [if)), one 
field with dimension (A) and 3 fields with dimension —2, (0, Xi 111 the case of 
perpendicular flow we have again if), B with dimension and x' with dimension —2. 
We see that the solution accommodates both cases because of the freedom we have 
in the definition of the physical fields. We define the physical fields (f> p h, ip p h as in the 
previous section (for simplicity we do not consider non-diagonal combinations at this 
point), and A ph = Ei' ~j V b v , B ph = Ei' h'h', Xph = Ei' Ci'0i'> <P P h = Ei' Vi'fa' and 
Xph = Ei' $i'0i'- The additional fusion rules involving by that we need are by x by = 
Efc' E$ f (f) k ,, b v x ipj = Efc' Dy^k' 

Previously, we discussed the constraint on the coefficients ati,f3j> originating from 
condition that the OPE of ip p h with itself give 4> p h as a minimal dimension field. This 
requirement is an integral part of the conformal approach to solving the Hopf equations. 
In the present case we have similar conditions which must be imposed on the coefficients 
occuring in the definitions of A p h, --Xph- Specifically we have 

A ph X A ph = Pph^Yl li'li' E i'f = Vk> 

i'j' 

Aph x ip ph = Xph->Y1 H' a 3 D ij = Cfc' 

i'j 

Bph x ip ph = Xph - = °k' (32) 

i'j 

So if we know on and 7^ we can find the rest of the coefficients in the case of the 
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parallel flow. The same holds for the perpendicular case if we know ctj and 5i>. In 
fact we expect that a change in these parameters is related to changing the boundary 
conditions that leave the spectrum invariant and so represent a kind of universality. 
These are the free parameters of the solution. 

At this stage we wish to make a few comments concerning solutions presented. It is 
apparent from Flohr's work [11] and the results presented here, that taking copies of a 
conformal model provides a way to incorporate the infrared condensate. If the copies 
are different then we have an asymmetry as can be seen from the following. Take as 
another possible solution c^i ® Cq^. Then we have two A = —2 and two A = — 1 fields, 
0i = *(-i|-i), 02 = ^i = *(-i|o)> ^2 = *(5|-i)- However there is only 

one A = field, namely b = It follows that b x b = 0i and so this naturally 

provides another solution of the case of perpendicular flow only. There may be other 
asymmetric tensor product solutions of this type. 

3 Conclusions 

In this paper we have considered the possibility of describing solutions of d = 2 ideal 
MHD turbulence within the framework of Flohr's logarithmic field theory solution in 
the pure fluid case. We have seen that the c 6j i <g> C6,i theory has sufficient degeneracy 
of dimension 0,-1 and —2 fields to be able to describe the additional fields one finds 
in the MHD case. Moreover the solutions of the corresponding Hopf equations was 
seen to be a consequence of the same hidden continuous symmetry that was exploited 
in the ordinary fluid case. This picture provides an interesting viewpoint of turbulent 
2-d fluids and plasmas as corresponding to different regions of parameter space within 
a single model. 

It would also be interesting to clarify further the exact relationship between the 
values of the real parameters entering the definitions of the physical fields ip p h, <f> p h, A p h 
and B ph given earlier, and boundary conditions that would be appropriate in a given 
physical situation. This question was already touched upon in [11]. Also one can 
imagine that our results should be extendible to the case of "dyonic" or duality invariant 
MHD [20 ], whose minimal model conformal solutions have been investigated in [21]. 

Finally, there has been some interest in the past, in studying perturbations of CFT 
solutions to turbulence which break the conformal symmetry [22-24]. The application 
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of such perturbations in LCFT's has been considered in [13] and [25] (and more recently 
in [27]). It would be interesting to apply these ideas to the particular MHD models we 
have considered here. 
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